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Abstract 

> 

We study the kinetics of the distribution function for charged particles of hard 
momentum in scalar QED. The goal is to understand the effects of infrared 



o 
in , 

CN | divergences associated with the exchange of quasistatic magnetic photons in 

the relaxation of the distribution function. We begin by obtaining a kinetic 

■ transport equation for the distribution function for hard charged scalars in a 

■ perturbative expansion that includes hard thermal loop resummation. Solv- 
Oh, ing this transport equation, the infrared divergences arising from absorption 
Q_i! and emission of soft quasi-static magnetic photons are manifest in logarithmic 

secular terms. We then implement the dynamical renormalization group re- 
summation of these secular terms in the relaxation time approximation. The 
distribution function (in the linearized regime) is found to approach equi- 
librium as 5nk{t) = Srik(t ) e ~' 2aT ( t - t o)M(i~'to)p] ; w ith /2 « uo p the plasma 
frequency and a = e 2 /4ir. This anomalous relaxation is recognized to be the 
square of the relaxation of the single particle propagator, providing a gener- 
alization of the usual relation between the damping and the interaction rate. 
The renormalization group approach to kinetics reveals clearly the time scale 
t re i ~ (aTlnfl/a]) -1 arising from infrared physics and hinges upon the sep- 
aration of scales t re i » ui^ 1 . 



I. INTRODUCTION 

The possibility of probing the quark-pluon plasma at the forthcoming RHIC at BNL and 
LHC at CERN, has sparked an intense activity in the understanding of collective excitations 
in ultrarelativistic plasmas and their potential experimental signatures (for reviews see 



1 



T0|). An important issue in this program is a consistent assessment of the time scales for 



thermalization and equilibration of the quasiparticles and their distribution functions, i.e. 
the damping and interaction rates [PHT3|j. At high temperatures usual perturbation theory- 



breaks down, but a consistent resummation program developed by Braaten and Pisarski [14 



|T8| provides a systematically improved perturbative expansion known as the hard thermal 



loop resummation (HTL). The HTL resummation incorporates screening corrections in a 
gauge invariant manner and is sufficient to render finite the damping rate of excitations at 
rest in the plasma [[HJ and transport cross sections [2(J . 



However these corrections are not sufficient to cure the infrared sensitivity of the damping 
rate of hard charged excitations which is dominated by the exchange of quasistatic magnetic 
photons and or gluons ||2lU22| . Whereas in QCD it is expected that a non-perturbative 
magnetic mass will provide an infrared cutoff and ameliorate the infrared sensitivity of the 
damping rate for hard charged particles |f2l|j22]| , in QED the transverse photons do not acquire 
a magnetic mass, screening is only dynamical through Landau damping and therefore the 
infrared singularities remain, possibly to all orders. This infrared sensitivity in the case 
of QED has led some authors to question the validity of the quasiparticle description and 
exponential relaxation that is the main concept behind the calculation of the damping rate 

ma 

These questions were recently clarified by the implementation of a Bloch-Nordsieck 
(eikonal) resummation of the infrared divergent diagrams in QED [^,^(| which leads to an 
anomalous real-time relaxation of the electron propagator of the form »s exp[— a T t ln(i u) p )]. 
This result has been recently confirmed in scalar QED (SQED) by implementing a dynamical 
renormalization group resummation directly in real time p^] . The advantage of this dynam- 
ical renormalization group approach is that it provides a real-time description of relaxation 
bypassing the limitation of the quasiparticle interpretation. In reference p7| the equivalence 
of the dynamical renormalization group to well studied situations at zero temperature was 
established: for cases with infrared threshold singularities at zero temperature it is equiva- 
lent to the Bloch-Nordsieck and the usual Euclidean renormalization group resummation. In 
the absence of infrared threshold singularities it leads to the real time evolution obtained by 
more conventional methods (for a more complete comparison the reader is referred to ||27|| ). 
For SQED at finite temperature it leads to the relaxation of the single particle propagator 
which agrees exactly with the results of the eikonal (Bloch-Nordsieck) approximation pH| in 
QED. This is one more example of the similarities between QED, QCD and SQED [fj8| in 
lowest order in the HTL resummation. 

In this article we explore a different but related question: since in hot QED and in SQED 
the infrared divergences associated with the emission and absorption of quasi-static magnetic 
photons lead to anomalous, non-exponential relaxation of the propagator of charged excita- 
tions p5|-f27|, what is the kinetic equation that describes the relaxation of the distribution 
function for these excitations?. This kinetic equation would be the equivalent of the Boltz- 
mann equation for the distribution function. 

When the quasiparticle picture is valid, the relaxation time approximation leads to the 
linearized equation 5Nk(t) = — 7^ 5Nk(t) with SNk(t) is the departure from equilibrium of 
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the quasiparticle distribution function and the relaxation rate 7& (inverse of the relaxation 
time) is time independent and simply related to the damping rate for single quasiparticles 
Tfc: Ik = p|-|TTf. If the single quasiparticle propagator is not exponentially damped it 
is reasonable to expect that the linearized kinetic equation for the distribution function will 
require a time dependent generalization of 7^. 

The goals of this article are a): a derivation of the proper kinetic equation for the case of 
hard charged (quasi) particles, b) interpretation of the relaxation time approximation and 
comparison to the relaxation of the single (quasi) particle Green's function. 

We study these issues within the context of SQED to make contact with previous results 
|2T| . This model is relevant to study the physics of the QGP because to lowest order in a 
SQED has the same infrared divergence and HTL structure as both QED and QCD [p8fl . 

The strategy: Our strategy leading to the kinetic equation follows closely the derivation 



of the Boltzmann equation presented in |2{J|3(J. It begins by defining a suitable number 



operator A 7 *.. In the case under consideration that of a hard scalar with momentum k > T 



the collective mode and the (renormalized) particle are indistinguishable pSyTOf , thus the 



number operator is the usual operator associated with asymptotic states. The time derivative 
of this operator is obtained from the Heisenberg equations of motion, and its average over an 
initial density matrix is performed using perturbation theory in terms of the non-equilibrium 
propagators I^^Oj. We find that the time evolution of iV& obtained in this perturbative 
expansion contains logarithmic secular terms in time. The dynamical renormalization group 
program is invoked to resum these secular terms. The resummed distribution function obeys 
a dynamical renormalization group equation which is recognized to be the generalization of 



the Boltzmann equation pi 



The results: There are two main results of our study, one of general scope and the 
other particular to the anomalous kinetics of hard charged particles in SQED (but likely 
to be shared by QED and QCD in lowest order in HTL resummation): a) we propose 
a microscopic derivation of quantum kinetic equations beginning from a non-equilibrium 
perturbative expansion in real time and using the dynamical renormalization group to resum 
the secular terms. This formulation allows to include other resummation schemes such as 
HTL directly in the derivation of the kinetic equation, b) More specifically to the problem 
of the kinetics of the distribution function for hard charged particles, we have focused on 
SQED as a model that bears many of the relevant features of QED and QCD. We find 
that the linearized version of the kinetic equation (relaxation time approximation) for this 
distribution function evolves in time as Snk(t) = 5rik(t ) e - 2aT (*-*o) in /*(*-*<>) f or ^ f g 
(with a = e 2 /(47r) the 'fine structure constant' and /2 ~ uo p ). This must be compared to 
the real-time evolution of the hard scalar propagator Gk{t) ~ e ~-aTtinpt 1 27 ] which reveals 
a similar relation between the relaxation time scales of the single particle Green's function 
and that of the distribution function as in the usual quasiparticle picture but in a manner 
that is not associated with pure exponential relaxation. 

A similar anomalous relaxation has been found by |25] in spinor QED via the Bloch- 
Nordsieck (eikonal) approximation. 

The article is organized as follows: in section II we introduce the theory and describe 
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the perturbative framework that leads to the kinetic equation to 0(a). The different con- 
tributions are interpreted as in-medium processes. In section III we focus on the relaxation 
time approximation (linearized kinetic equation), recognize the secular terms arising in the 
perturbative expansion and implement the dynamical renormalization group resummation 
of these secular terms in the asymptotic long time limit. In section IV we provide an inter- 
pretation of the resummation scheme resulting from the dynamical renormalization group 
and generalize the solution to include short time transients. Section V summarizes our con- 
clusions, provides an assessment of the potential impact of these results and poses other 
questions. 



II. THE KINETIC EQUATION 

We propose to study the relaxation of hard charged scalars in scalar QED as a prelude 
to studying the more technically involved cases of QED and QCD. Scalar QED shares many 
of the important features of QED and QCD in leading order in the HTL resummation 
|28fl . Furthermore, the infrared physics in QED captured in the eikonal approximation 
(Bloch-Nordsieck) as clearly explained in |25[ has been reproduced recently via the dynamical 



renormalization group in SQED [p7|| , thus lending more support to the similarities of both 
theories at least in leading HTL order. 

In the Abelian theory under consideration, it is rather straightforward to implement a 
gauge invariant formulation by projecting the Hilbert space on states annihilated by Gauss' 
law. Gauge invariant operators can be constructed and the Hamiltonian and Lagrangian 
can be written in terms of these. The resulting Lagrangian is exactly the same as that in 
Coulomb gauge |30| (for more details see p2| ) and is given by 

£ = «9 M $t _ m 2 $ t$ + -d^l T ■ d»A T - el T ■ j T - e 2 A T ■ A T $ f $ + 
+ l -(VAvf + e 2 Al $t$ + eA oP , 



j T = i($ty T $ - Vr* 1 "*) ; p = -i ($$ f - $ f $) 



where we have traded the instantaneous Coulomb interaction for a gauge invariant Lagrange 
multiplier field Aq which should not be confused with a time component of the gauge field. 
A? is the transverse component satisfying V ■ Aj-(x,t) = 0. Since we are only interested in 
obtaining the infrared behavior arising from finite temperature effects we do not introduce 
the renormalization counterterms to facilitate the study, although these can be systematically 
included in our formulation [|27[| . The finite temperature behavior is ultraviolet finite. The 
non-equilibrium generating functional requires the fields on the forward and backward time 
branches as discussed in |30| (for more details on the non-equilibrium formulation see |33 



We anticipate that the Coulomb interaction will not be relevant to the infrared physics 
because the longitudinal photon will be screened with a screening mass m s oc eT. Only 
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the transverse photons will lead to infrared divergences pi|,p5|,pT| an d therefore we neglect 
the contribution from longitudinal photons (the Coulomb interaction). Furthermore we will 
consider a neutral system with vanishing chemical potential. It is convenient to perform a 
spatial Fourier transform at some given time and write 



d 3 k 
(2tt) 3 / 2 
d 3 k 



4>{k) e 



ik-x 



4>{k) e 



ik-x 



i 



V^k 

2 



(2vr) 3 / 2 

And the spatial Fourier transform for the transverse gauge field 

d 3 k 



ak 



A T x) 



(2tt) 3 / 2 



A T (k) e 



\k-x 



Since we are interested in hard charged massive scalars, uo k ~ k the collective modes are indis- 
tinguishable from the single particle excitations [25, II]]. The soft collective modes require a 
definition of the quasiparticle number operator, for which the formulation recently proposed 
in terms of quasiparticles PB| could prove useful. This case will be explored elsewhere. 

The number of positively charged scalars (which at zero chemical potential is equal to 
the number of negatively charged scalars) is then given by the following Heisenberg operator 



n + (k,t) = a\(t)a k (t) 



2uJk 



ft(k,t)j>(k,t)+u 2 k ^{k,t)(f>{k,t) 



+ iuo k t)<fi(k, t) - f (fc, t)(f)(k, t) 

Using the Heisenberg equations of motion we obtain 

d . \ , + J f d 3 q 



h + (k,t) 



e 



e 

+ — 

UJk 



'd_ 
dt 



(2tt) 3 / 2 



iu k j (j)(k,t) 
K T {q) ■ A T (q,t)<f){k - q,t) 



K T (q) =k-q(k-q) 



We need to take the expectation value of this Heisenberg operator in the initial density 
matrix. The following identities prove useful 



Tr [A{t)B{t)C{t)p{U)} = Tr [B{t)C{t) P {U)A{t)} 

= Tr [C(t)p(U)A(t)B(t)} 



(B+(t)C+(t)A-(t)) 
(C+(t)A-(t)B-(t)) 



where the ± correspond to operators defined on the forward (+) and backward (— ) time 
branches, which are obtained as functional derivatives with respect to sources of the forward 
and backward time evolution operators |33|]- [ 35| . These identities allow to extract the time 
derivatives outside of the expectation values and to avoid the potential Schwinger terms 
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associated with time ordering. Now the expectation value of n + (k,t) given above can be 
written as follows 



(h+(k,t)) 



7 - Wk 



e_ (d_ 

LU k \dt 

u k \dt' lUk / 



d 3 q 



(2tt) 3 / 2 
(2vr) 3 /2 



%(?)(^T(?,^ tr (H9,^ + (M')> \t=V + 
K ilT (Q) (Avr(£ t)<P + (k - q, t)^'~ (k, t')> \ t=v 



This expectation value is computed in non-equilibrium perturbation theory in terms of the 
non-equilibrium propagators and vertices. Since we are interested in the infrared region in 
the internal loop momenta, we must consider the HTL resummed photon propagator, but 
the scalar propagator need not be resummed because k > T and the scalar in the loop is 
hard and massive. The required non-equilibrium propagators are the following |30|,27,|35| 
•Scalar Propagators (zero chemical potential) 



where (a, b) e {+,—}. 



G++(t, = G>(t, t')B(t - 1') + G<(t, t')B(t' - 1) 
G p -(t,t') = G>(t,l/)e(t' -t) + G<(t,t')Q(t - 
G^(t,t')=G<>\t,t'), 

G>(t,t') 



2uo r 



2uo r 



n p e 



-iujp(t—t') 



+ (1 + n v ) e 



iw p (t—t') 



n, 



»HTL dressed Photon Propagators |]T0] , pT 



(A$(q,t)A$](-q,t')) = -iG%(q;t,t f ) 

g± + (q; t, t') = ViM [g q >(t, t')e(t - ?) + g<(t, t')Q(t - 1) 

gr j -(q;t,t')=V lJ (q) [0\ 'MW " *) + 0? (*,0©(* " - 

g^(q;t,t') = Va%g<^ m, 

g>(t, = ~ / d 1o PT(q ; q) [(1 + Nje-**>W + 
g<(t, t')= l -( dq p T (q ; q) K + (1 + NJ 



N. 



<lo 
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Where we have used the properties ||2"7|,|TC| 

Pr(-Qo,q) = -Pt(Qo,q) 
The HTL spectral density is given by [p7| , 

Pr(qo,q) 

S/(g ,g) 
^ (go, q) 



1 


£/ 


71 [ql - 


-<? 2 - 


ne 2 T 2 


9o / 


12 


? V 


e 2 T 2 


9 ?0 


12 


g z 



jV(-g ) = -[l + iV(g )] 



+ sign(go) <%o ~ 



1-^ 

<7 2 , 



(2.1) 

(2.2) 
(2.3) 



+ 



q 



9o - 9 



where a> p (g) is the plasmon pole and Z(q) its (momentum dependent) residue, which will 
not be relevant for the following discussion. The important feature of this HTL resummed 
spectral density is its support below the light cone, i.e. for q 2 > q 2 ,, the imaginary part 



originates in the process of Landau damping |2l] from scattering of quanta in the medium. 
Above Vij(q) is the transverse projection operator: 



<Mi 
q 2 



Using the above expressions for the non-equilibrium propagators, and after some tedious but 
straightforward algebra, we find the expectation value (h k (t)) to lowest order in perturbation 
theory 0(e 2 ) is given by 



("*(*)> 



e 2 k 2 



4:7T 2 UJk JO iOt, 



q 2 dq 



k+q 



dcosO (1 — cos 2 



-i 



{ [(1 + N Qo )(l + n %+ -){l + n k ) - N qo 
N a (l + n 



dq a pT{q ,q) [ dt' x 
Ju 

cos(g Q + LU£ + - + uJk)(t - t') + 



+ n h ) - (1 + N qo ) n Uq ,n k \ cos(-g + uJ^+u k )(t - t') + 
N io (1 + n k ) - (1 + NJ(1 + n u? ) n k ] cos(g + Uj^-u k )(t - t')} (2.4) 

where 9 is the angle between k and q. The different contributions have a very natural 
interpretation in terms of gain — loss processes. The first term in brackets corresponds to 
the process — > 7* + s + s minus the process 7* + s + s — * 0, the second term corresponds to 
7* — > s + s minus s + s — > 7*, and the last term corresponds to the scattering in the medium 
7* + s — > s minus the inverse process s —> 7* + s where 7* refers to the HTL- dressed photon 
and s, s refer to the charged quanta of the scalar field $. 



III. RELAXATION TIME APPROXIMATION: SECULAR TERMS 

We will assume that there is an equilibrium solution for the distribution function and 
that at time t = to the distribution function for a fixed mode k is disturbed slightly off 
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equilibrium so that n k (t < t Q ) = n e k g ; n k {t = t ) 



n 



eq 



+ Sn k (to) while the rest of the modes 



remain in equilibrium. We want to study the time evolution of the perturbed distribution 
in the linearized approximation. This leads to a kinetic equation ( |2.4| ) that is linear in Sn k , 
since consistent with the condition that only the mode of wavevector k is off equilibrium, 
we set Sn^ + ^ = for q ^ and the point q = does not contribute because the integration 
measure vanishes there. 

This is known as the relaxation time approximation. Only in this approximation the 
relaxation of the distribution function for the number of particles is related to the relaxation 
of single particle Green's function. Since the propagators entering in the perturbative expan- 
sion of the kinetic equation are in terms of the distribution functions at the initial time the 
time integration can be done straightforwardly leading to the following linearized equation 

c 2 k 2 q 2 dq 

ShJt) = 5n k (t ) — — / / (1 - cos 2 6>) dcos6 / dq Q p T (g Q ,g) x 

All UJh JO U}r,-fJ-l J-oo 



l + N qo + n 



k+q 

sin(g, 



to) 



l + N qo + n 



b k+q 



q + u 



k+q 



(3.1) 



where we have used the properties ( (2.1| ) to combine the first and second terms in the kinetic 
equation (|2.4|). 

Before proceeding further, it is illustrative to analyze the above equation in the long time 
limit t » t Q . If the denominators in ( |3.1| ) have isolated zeroes in the region of integration we 
can use the approximation sin[f2r]/f2 t ps°° ttS(Q) which is the usual approximation leading 
to Fermi's Golden rule. Thus if 'there are no singularities arising from the integrals and if the 
resulting delta functions are satisfied one finds a linear secular term in time upon integrating 
in time the rate equation (|3.1|) . This is a perturbative signal of pure exponential relaxation 
at long times |HJ. However in the case under consideration there are infrared singularities 



in the integrand p5|j27|1 and the long-time limit must be studied carefully. 
For this purpose it is convenient to introduce the following spectral density 



z 2 k 2 



Ti 



q 2 dq 



LOt 



'1 — cos 9) dcos9 



-l 



h pT(q Q , q) x 



1 + N(q ) + n n+? ] 5(uj-q, 



to 



k+q> 



(3.2) 



which is the same as that studied within the context of the relaxation of the amplitude of a 
mean field in E7[ and in the eikonal approximation p5|. 



In terms of this spectral density we obtain the time derivative of distribution function in 
the form 



5h k (t) = -a T k (t) 5n k (t ) 
1 



duj p(k; uj) 



sin [(a; - u k )(t - t )} 
(to - u k ) 



+ (uj 



-UJ 



a 



e 

47T 



(3.3) 
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which upon integrating in time with initial condition at t a leads to the form 



5n k {t) = 5n k {t ) - a J* T k {t') dt' 



We focus on the case of hard external momemtum k pa T but on the infrared region of the 
loop integral q Q , q < eT in the spectral function ( |3.2| ) f2"5|]2"7|| . This is the region dominated by 
the exchange of soft, HTL resummed transverse photons [0,0 and that dominates the long 
time evolution of the distribution function. In this region we can replace uj k ~ k ; uJ k+ ^ ~ 
k + q cos 6. 

Potential secular terms (growing in time) could arise in the long time limit t » t 
whenever the denominators in ( |3.3|) vanish, i.e. for the region of frequencies uj ~ ~ k and 
uj pa —u k pa —k. For uj pa uj k pa k we see that the argument of the delta function in (3.2) 
vanishes in the region of the Landau damping cut of the exchanged transverse photon q 2 < q 2 
and contributes to the infrared behavior. On the other hand, for uj pa —uj k ~ —k the delta 
function in ( |3.2| ) is satisfied for q Q pa —2k, and this region gives a negligible contribution to 
the long time dynamics. Therefore only the first term in ( |3.3|) (with uj — uj k ) contributes in 



the long time limit. This term is dominated by the Landau damping region of the spectral 
density of the exchanged soft photon given by ( |2.2| ) since for uj ~ uj k the argument of the 
delta function is q Q + q cos(6 l ) and this is the region where the imaginary part (12.3] ) has 
support. The second contribution (with uj + uj k ) oscillates in time and is always bound and 
perturbatively small. 

In the hard limit u k pa k, and for q \ q < eT << T we neglect the contribution from n<£ +( ~- 
(hard and massive scalar) and replace N(q ) ~ T/q a . The spectral density is found to be 
given by |2"7 



p{k;uj) 



♦ fc e 2 kT 



In 



7T" 



uj — k 



+ 0{uj - k) 



where /i is an infrared cutoff /x 

In the limit fi(t — t ) >> 1 we find [[27 



^ uj p pa 



eT [0,0,0. 



to 



T k (t')dt' 



n(t-t )»i 



2T (t — t Q ) ln[/i(t — t Q )] + non-secular terms 



(3.4) 



with pL = fiexpl'jE — 1] and is the Euler-Mascheroni constant. In lowest order in per- 
turbation theory, the occupation number that enters in the loops are those at the initial 
time. Obviously perturbation theory breaks down at time scales t — t Q pa |aTlna| _1 where 



we used fx 



0J r . 



eT. This situation is similar to the case of a weakly damped harmonic 



oscillator when a perturbative solution in terms of the damping coefficient is sought. Such 
a perturbative solution has secular terms that grow linearly in time in lowest order, which 
reflect the perturbative expansion of the damping exponential, i.e. e~ 7 * ~ 1 — yt -\ — -. This 
perturbative expansion breaks down at time scales 0(1/7). The dynamical renormalization 
group is a systematic generalization of multi-time scale analysis and sums the secular terms, 
thus improving the perturbative expansion |J7| . For a discussion of the dynamical renormal- 
ization group in other contexts, including applications to quantum field theory problems see 
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We now implement the dynamical renormalization group to sum the secular divergences 
and to improve the perturbative expansion using the formulation advanced in |27j . 

To achieve this purpose we introduce a renormalization constant for the occupation num- 
ber that absorbs the secular divergences at a fixed time scale r and write 

5n k (t ) = 5n k (r) Z(r,t ) ; Z(t, t Q ) = 1 + a z x (t, t Q ) H (3.5) 

and request that the coefficients z n cancel the secular divergences proportional to a n at a 
given time scale r. To lowest order the choice 

zi(r,t ) = f T T k (t') dt' (3.6) 

leads to the renormalized distribution function at time t in terms of the updated occupation 
number at the time scale r 

8n k {t) = 6n k {r){l-a f ' T k (t') dt' 



However, the occupation number Sn k (t) cannot depend on the arbitrary renormalization 
scale r, this independence on the renormalization scale leads to the renormalization group 
equation to lowest order: 

^^- + aT k (T)5n k (r) = 0. (3.7) 

This renormalization group equation is now clearly of the form of a Boltzmann equation in 
the relaxation time approximation with a time dependent rate. 

Now choosing the renormalization scale to coincide with the time t in the solution of 
( p.7|) as is usually done in the scaling analysis of the solutions to the renormalization group 
equations, we find that the distribution function in the linearized approximation evolves in 
time in the following manner 

5n k (t) = 5n k (t )e- a Si Tk{t ' )dt ' (3.8) 

with the initial conditions 

6n k (t = t Q ) = 5n k (t ) ; 5h k (t) \ t=to = (3.9) 

in agreement with the perturbative expression fl3.1| ). In the long time limit ~p(t — t Q ) » 1 



and using ( |3.4j) we find that the distribution function relaxes towards equilibrium as 

Sn k (t) M "= )>>1 6n k (t )e-* oa '<- t - t °W t - t °M (3.10) 



The exponent in (|3.10|) must be compared to that found for the relaxation of the mean field 
(&(k,t)) in 1^3- We thus find that the distribution function approaches equilibrium as the 
square of the mean field, a result which is in agreement with the usual relation between the 
damping and the interaction rate. 

Furthermore, the dynamical renormalization group resummation of the perturbative ex- 
pansion reveals a relaxation time scale for the distribution function of the hard charged 
particles given by 

trei ~ laTlnaT 1 . 



10 



IV. INTERPRETATION OF THE RG IN KINETICS 



The interpretation of eq. (|3.4j ) and the renormalization of the distribution function given 
by eq.( p.5| ) are clear: having prepared the initial distribution at a time t a we evolve the 
distribution forward in time using perturbation theory but during a time scale r such that 
the perturbative expansion is still valid, i.e. t Td » (r—t ). Secular terms begin to dominate 
the perturbative expansion at a time scale (r — t ) » ujp 1 , if there is a separation of time 
scales such that t re i » (r — t a ) » perturbation theory is reliable in this regime but 
secular terms appear and can be isolated. Thus the perturbative expansion is valid for a 
large time scale, at a given renormalization scale r within this interval the result of the 
perturbative expansion is to 'reset' the occupation number, from 5nj.(to) to Sn^r). Having 
'reset' the occupation number at this time scale in which perturbation theory is reliable, 
this new occupation number is used as an initial condition to iterate forward in time using 
perturbation theory to another time that again is within the perturbatively reliable region. 

The renormalization condition ( |3.5[ ) is precisely the resetting of the occupation number. 
Now the perturbative expansion has been improved and can be extended to longer time scales 
by using the updated occupation number as an initial condition at the scale r. The dynamical 
renormalization group equation ( |3.7| ) is the differential form of the procedure of 'resetting' 
followed by perturbative evolution in real time. This renormalization group equation (|3.7| ) 



is identified with the Boltzmann equation in the relaxation time approximation |pl |. 

We note that there is certain freedom in the choice of the 'renormalization counterterms' 
z n ( p.5| ): rather than choosing Z\ as in eq.( |3.6| ) we could have chosen to cancel only the 
secular part given by ( |3.4|) rather than the whole integral. In the language of renormalization 
this choice would correspond to choosing the counterterms to cancel only the divergent 
contribution, i.e. we could have chosen 

z x = 2T(r - Q ln[/x(r - t )\ (4.1) 

which would lead to the renormalization group equation 

^J^l + 2aT {ln[//(r - t )] + 1} 8n k (r) = (4.2) 

Now choosing the renormalization scale r with the time t in the solution of ( f4.2| ) we 
would find the asymptotic form 

5n k (t) = 5n k {t ) e -2aT(t-t o) M(*-*o)A] 

Obviously the two choices ( |3.6| ) and ( ^4.1| ) lead to the same asymptotic behavior, because 
they only differ by non-secular terms which are subleading in the asymptotic regime. The 
choice ( |3.6| ) also contains non-secular contributions which are bound in time and therefore 
always perturbative, i.e. finite terms in the usual renormalization program. The advantage 
of choosing fl3.6|) is that the solution (|3~8| ) includes transient effects and obeys the initial 
conditions ( |3.9|) explicitly. 
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In this particular case, the validity of the renormalization group approach hinges upon 
a separation of time scales: the time scale at which the occupation number is reset, i.e. r 
has to be close to the scale t so that the perturbative expansion has been improved but not 
too long so that the occupation number has changed much during this scale. Since the scale 
of the kernel is given by the plasma frequency uo p pa eT which enters in the self-energy of 
the HTL resummed photon and the infrared cutoff /i, the separation of scales is guaranteed 
in weak coupling since t re i ~ |aTln(a)p 1 >> co^ 1 ~ (eT) -1 . This is a general statement 
for the validity of any kinetic approach, that is, there must be a clear separation of scales 
between the relaxation and the microscopic time scales. 

This interpretation of kinetics and the Boltzmann equation via the dynamical renor- 
malization group provides a natural answer to the question: What does the Boltzmann 
equation sum??. In a typical kinetic approach one writes a Boltzmann equation by com- 
puting in perturbation theory some transition matrix element which enters in an integral 
kernel. Obviously the full solution of the Boltzmann equation is non-perturbative even when 
the transition matrix element has been computed in perturbation theory. In particular in 
the relaxation time approximation the coupling constant appears in the solution in an ex- 
ponential form via the relaxation time or damping rate. What our analysis in terms of the 
renormalization group reveals is that what is being summed by the Boltzmann equations are 
secular terms in time. 

In this article we have focused on the derivation of the kinetic equation and the renormal- 
ization group resummation in the relaxation time approximation to compare the relaxation 
of the distribution function of hard charged scalars to that of the single particle Green's 
function [ERE71]. We relegate the treatment of the full non-linear kinetic equations via the 



dynamical renormalization in its full generality to a forthcoming article [31 



V. CONCLUSIONS, IMPLICATIONS AND MORE QUESTIONS: 



Our goal with this article was to understand the relaxation of the distribution function 
of charged fields with hard momenta in a hot gauge theory. The detailed study of the relax- 
ation of the single charged (quasi)particle Green's function in the eikonal (Bloch-Nordsieck) 
approximation in references [^,[2(J revealed an anomalous relaxation as a consequence of 
the infrared divergences associated with the emission and absorption of HTL resummed 
transverse photons. This novel form of relaxation of the single particle Green's function was 



confirmed in terms of a dynamical renormalization group resummation in |p7 
There are several important implications of these results: 



Although in QED there is no magnetic mass, the analisis of |p5|-|27| reveals that the 
time scale for relaxation of fast moving single particle charged excitations is given by 
2|aTln(a)|~ 1 . This is the same time scale that arises in QCD after resummation and 
assuming that a magnetic mass m mag « rriD ~ gT cuts off the infrared divergences 
in the damping rate | 2~I| , |2~2 |. Thus the conjecture in [[H] is confirmed through a non- 



exponential relaxation. The advantage of a real time analisis is that the time variable 
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serves as an infrared cutoff. Furthermore the renormalization group approach is not 
restricted to the quasiparticle picture. 



• The distribution function of fast charged excitations in the linearized approximation 
relaxes also with an an anomalous, non-exponential form that asymptotically is the 
square of the single particle Green's function. Therefore this provides another con- 
firmation of the relation between the 'damping' and 'interaction' rates but with non- 
exponential relaxation. The relaxation time scale for the distribution function is given 
by t re i « |aTln(a)| _1 , i.e. one half of the relaxation time scale for the single particle 
Green's function. 

• This result suggests that at very large energy when the QCD coupling is small, i.e. in 
the initial stages of an Ultrarelativistic Heavy Ion Collision, the thermalization rate 
of quarks could be as fast as that of gluons. The in-medium effects that lead to the 
logarithmic enhancement of the thermalization scale could compensate for the group 
factors that result in a larger thermalization rate for gluons. 

There are still unresolved and important questions that deserve attention. In particular 
the treatment of collective modes. Here we have focused on hard excitations, and in the limit 
of large momenta k » eT collective and single particle excitations are indistinguishable. 
Hence the definition of the number of charged excitations coincides with that for asymptotic 
single particle states. However for k < eT we must treat the distribution function of the 
collective modes in terms of an interpolating operator that counts these collective modes. 
Furthermore for soft collective modes, not only the internal photon but the charged particle 
propagator and the vertices must be HTL resummed 0,^] and this situation must be 
studied in detail. We hope to report on progress on some of these issues along with a 
treatment of non-linear kinetics and resolution of pinch singularities with the dynamical 



renormalization group in a forthcoming article [31 
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